ebab7d0da53c51438720d60c6e2e467c - /page/wikibook/entry.mu

<!DOCTYPE html>

<title>Standard step method</title>

</head>

<span id="openzim-page-title" class="mw-page-title-main"><span class="mw-page-title-main">Standard step method</span></span>

</h1>

</header>

<p>The <b>standard step method (STM)</b> is a computational technique utilized to estimate one-dimensional surface water profiles in open channels with gradually varied flow under steady state conditions. It uses a combination of the energy, momentum, and continuity equations to determine water depth with a given a friction slope <span class="mwe-math-element mwe-math-element-inline"><span class="mwe-math-mathml-inline mwe-math-mathml-a11y" style="display: none;"><math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle (S_{f})}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle (S_{f})}</annotation>

</semantics>

</math></span><img src="./fb53cb3f3e0bc326021492f4d366dff2e27105d1.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:4.371ex; height:3.009ex;" alt="{\displaystyle (S_{f})}" loading="lazy"></span>, channel slope <span class="mwe-math-element mwe-math-element-inline"><span class="mwe-math-mathml-inline mwe-math-mathml-a11y" style="display: none;"><math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle (S_{0})}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle (S_{0})}</annotation>

</semantics>

</math></span><img src="./6f5bd4ceefda1801d074f796dade72b59cca9035.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:4.288ex; height:2.843ex;" alt="{\displaystyle (S_{0})}" loading="lazy"></span>, channel geometry, and also a given flow rate. In practice, this technique is widely used through the computer program <a href="HEC-RAS" title="HEC-RAS">HEC-RAS</a>, developed by the US Army Corps of Engineers Hydrologic Engineering Center (HEC).<sup id="cite_ref-1" class="reference"><a href="#cite_note-1"><span class="cite-bracket">[</span>1<span class="cite-bracket">]</span></a></sup>

</p>

<div class="mw-heading mw-heading2"><h2 id="Open_channel_flow_fundamentals">Open channel flow fundamentals</h2></div>

<p>The energy equation used for <a href="Open_channel_flow" class="mw-redirect" title="Open channel flow">open channel flow</a> computations is a simplification of the Bernoulli Equation (See <a href="Bernoulli_Principle" class="mw-redirect" title="Bernoulli Principle">Bernoulli Principle</a>), which takes into account pressure head, elevation head, and velocity head. (Note, energy and head are synonymous in Fluid Dynamics. See <a href="Pressure_head" title="Pressure head">Pressure head</a> for more details.) In open channels, it is assumed that changes in atmospheric pressure are negligible, therefore the “pressure head” term used in Bernoulli’s Equation is eliminated. The resulting energy equation is shown below:

</p>

<mfrac>

<msup>

</mrow>

</msup>

<mrow>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle H=z+y+{\frac {v^{2}}{2g}}}</annotation>

</semantics>

</math></span><img src="./0bb59751df63b21088b4f2042cb0d12cb391bcbb.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.338ex; width:16.201ex; height:6.176ex;" alt="{\displaystyle H=z+y+{\frac {v^{2}}{2g}}}" loading="lazy"></span>            <big><b>Equation 1</b></big></dd></dl>

<p>For a given flow rate and channel geometry, there is a relationship between flow depth and total energy. This is illustrated below in the plot of energy vs. flow depth, widely known as an E-y diagram. In this plot, the depth where the minimum energy occurs is known as the critical depth. Consequently, this depth corresponds to a <a href="Froude_Number" class="mw-redirect" title="Froude Number">Froude Number</a> <span class="mwe-math-element mwe-math-element-inline"><span class="mwe-math-mathml-inline mwe-math-mathml-a11y" style="display: none;"><math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle (F_{n})}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle (F_{n})}</annotation>

</semantics>

</math></span><img src="./e28830549cbbd029f90ee9a4358c62c72455cead.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:4.522ex; height:2.843ex;" alt="{\displaystyle (F_{n})}" loading="lazy"></span> of 1. Depths greater than critical depth are considered “subcritical” and have a Froude Number less than 1, while depths less than critical depth are considered supercritical and have Froude Numbers greater than 1.

</p>

<msub>

</mrow>

</msub>

<mfrac>

<mrow>

<mfrac>

</mfrac>

</mrow>

<msup>

</mrow>

</msup>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle F_{n}={\frac {v}{(g{\frac {A}{B}})^{0.5}}}}</annotation>

</semantics>

</math></span><img src="./ddc3bd5589a6798e404e79e43ad0a290ad4090f8.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -3.505ex; width:13.99ex; height:6.343ex;" alt="{\displaystyle F_{n}={\frac {v}{(g{\frac {A}{B}})^{0.5}}}}" loading="lazy"></span>           <big><b>Equation 2</b></big></dd></dl>

<p>Under steady state flow conditions (e.g. no flood wave), open channel flow can be subdivided into three types of flow: uniform flow, gradually varying flow, and rapidly varying flow. Uniform flow describes a situation where flow depth does not change with distance along the channel. This can only occur in a smooth channel that does not experience any changes in flow, channel geometry, roughness or channel slope. During uniform flow, the flow depth is known as normal depth (yn). This depth is analogous to the terminal velocity of an object in free fall, where gravity and frictional forces are in balance (Moglen, 2013).<sup id="cite_ref-3" class="reference"><a href="#cite_note-3"><span class="cite-bracket">[</span>3<span class="cite-bracket">]</span></a></sup> Typically, this depth is calculated using the <a href="Manning_formula" title="Manning formula">Manning formula</a>. Gradually varied flow occurs when the change in flow depth per change in flow distance is very small. In this case, hydrostatic relationships developed for uniform flow still apply. Examples of this include the backwater behind an in-stream structure (e.g. dam, sluice gate, weir, etc.), when there is a constriction in the channel, and when there is a minor change in channel slope. Rapidly varied flow occurs when the change in flow depth per change in flow distance is significant. In this case, hydrostatics relationships are not appropriate for analytical solutions, and continuity of momentum must be employed. Examples of this include large changes in slope like a spillway, abrupt constriction/expansion of flow, or a hydraulic jump.

</p>

<div class="mw-heading mw-heading2"><h2 id="Water_surface_profiles_(gradually_varied_flow)">Water surface profiles (gradually varied flow)</h2></div>

<p>Typically, the STM is used to develop “surface water profiles,” or longitudinal representations of channel depth, for channels experiencing gradually varied flow. These transitions can be classified based on reach condition (mild or steep), and also the type of transition being made. Mild reaches occur where normal depth is subcritical (yn > yc) while steep reaches occur where normal depth is supercritical (yn<yc). The transitions are classified by zone. (See figure 3.)

</p><p><span class="mw-default-size" typeof="mw:File"></span>

</p><p><b>Figure 3.</b> This figure illustrates the different classes of surface water profiles experienced in steep and mild reaches during gradually varied flow conditions.<sup id="cite_ref-4" class="reference"><a href="#cite_note-4"><span class="cite-bracket">[</span>4<span class="cite-bracket">]</span></a></sup> Note: The Steep Reach column should be labeled "Steep Reach (yn<yc).

</p><p>The above surface water profiles are based on the governing equation for gradually varied flow (seen below)

</p>

<mfrac>

<mrow>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

</mrow>

</mrow>

</msubsup>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {dy}{dx}}={\frac {S_{0}-S_{f}}{1-F_{n}^{2}}}}</annotation>

</semantics>

</math></span><img src="./2a9340d7f3a14d91ee39e016d045aef84bccceda.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.505ex; width:15.197ex; height:6.343ex;" alt="{\displaystyle {\frac {dy}{dx}}={\frac {S_{0}-S_{f}}{1-F_{n}^{2}}}}" loading="lazy"></span>           <big><b>Equation 3</b></big></dd></dl>

<p>This equation (and associated surface water profiles) is based on the following assumptions:

</p>

<dl><dd><ul><li>The slope is relatively small</li>

<li>Channel cross-section is known at stations of interest</li>

<li>There is a hydrostatic pressure distribution</li></ul></dd></dl>

<div class="mw-heading mw-heading2"><h2 id="Standard_step_method_calculation">Standard step method calculation</h2></div>

<p>The STM numerically solves equation 3 through an iterative process. This can be done using the bisection or Newton-Raphson Method, and is essentially solving for total head at a specified location using equations 4 and 5 by varying depth at the specified location.<sup id="cite_ref-5" class="reference"><a href="#cite_note-5"><span class="cite-bracket">[</span>5<span class="cite-bracket">]</span></a></sup>

</p>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle H_{2}=H_{1}-h_{f}}</annotation>

</semantics>

</math></span><img src="./587ca8002731882d60b5b5f87163d3f39c183fcd.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:14.385ex; height:2.843ex;" alt="{\displaystyle H_{2}=H_{1}-h_{f}}" loading="lazy"></span>           <big><b>Equation 4</b></big></dd></dl>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle H_{2}=h_{vel}+h_{ele}}</annotation>

</semantics>

</math></span><img src="./269a99e0641a4bcf3ac211fb08865ce733f890e9.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:16.143ex; height:2.509ex;" alt="{\displaystyle H_{2}=h_{vel}+h_{ele}}" loading="lazy"></span>           <big><b>Equation 5</b></big></dd></dl>

<p>In order to use this technique, it is important to note you must have some understanding of the system you are modeling. For each gradually varied flow transition, you must know both boundary conditions and you must also calculate length of that transition. (e.g. For an M1 Profile, you must find the rise at the downstream boundary condition, the normal depth at the upstream boundary condition, and also the length of the transition.) To find the length of the gradually varied flow transitions, iterate the “step length”, instead of height, at the boundary condition height until equations 4 and 5 agree. (e.g. For an M1 Profile, position 1 would be the downstream condition and you would solve for position two where the height is equal to normal depth.)

</p>

<div class="mw-heading mw-heading3"><h3 id="Newton–Raphson_numerical_method">Newton–Raphson numerical method</h3></div>

</p><p>Computer programs like excel contain iteration or goal seek functions that can automatically calculate the actual depth instead of manual iteration.

</p>

<div class="mw-heading mw-heading3"><h3 id="Conceptual_surface_water_profiles_(sluice_gate)">Conceptual surface water profiles (sluice gate)</h3></div>

<p>Figure 4 illustrates the different surface water profiles associated with a sluice gate on a mild reach (top) and a steep reach (bottom). Note, the sluice gate induces a choke in the system, causing a “backwater” profile just upstream of the gate. In the mild reach, the <a href="Hydraulic_jump" title="Hydraulic jump">hydraulic jump</a> occurs downstream of the gate, but in the steep reach, the hydraulic jump occurs upstream of the gate. It is important to note that the gradually varied flow equations and associated numerical methods (including the standard step method) cannot accurately model the dynamics of a hydraulic jump.<sup id="cite_ref-6" class="reference"><a href="#cite_note-6"><span class="cite-bracket">[</span>6<span class="cite-bracket">]</span></a></sup> See the <a href="Hydraulic_jumps_in_rectangular_channels" title="Hydraulic jumps in rectangular channels">Hydraulic jumps in rectangular channels</a> page for more information. Below, an example problem will use conceptual models to build a surface water profile using the STM.

</p>

<div class="mw-heading mw-heading2"><h2 id="Example_problem">Example problem</h2></div>

</p><p><big><b>Solution</b></big>

</p><p><span typeof="mw:File"></span>

</p><p><br>

</p><p><br>

</p><p><br>

Using Figure 3 and knowledge of the upstream and downstream conditions and the depth values on either side of the gate, a general estimate of the profiles upstream and downstream of the gate can be generated. Upstream, the water surface must rise from a normal depth of 0.97 m to 9.21 m at the gate. The only way to do this on a mild reach is to follow an M1 profile. The same logic applies downstream to determine that the water surface follows an M3 profile from the gate until the depth reaches the <a href="Conjugate_depth" class="mw-redirect" title="Conjugate depth">conjugate depth</a> of the normal depth at which point a hydraulic jump forms to raise the water surface to the normal depth.

</p><p><b>Step 4:</b> Use the Newton Raphson Method to solve the M1 and M3 surface water profiles. The upstream and downstream portions must be modeled separately with an initial depth of 9.21 m for the upstream portion, and 0.15 m for the downstream portion. The downstream depth should only be modeled until it reaches the conjugate depth of the normal depth, at which point a hydraulic jump will form. The solution presented explains how to solve the problem in a spreadsheet, showing the calculations column by column. Within Excel, the goal seek function can be used to set column 15 to 0 by changing the depth estimate in column 2 instead of iterating manually.

</p><p><br>

</p><p><span typeof="mw:File"></span>

</p><p>Table 1: Spreadsheet of Newton Raphson Method of downstream water surface elevation calculations

</p><p><br>

<b>Step 5:</b> Combine the results from the different profiles and display.

</p><p><span typeof="mw:File"></span>

</p><p>Normal depth was achieved at approximately 2,200 meters upstream of the gate.

</p><p><b>Step 6:</b> Solve the problem in the HEC-RAS Modeling Environment:

</p><p>It is beyond the scope of this Wikipedia Page to explain the intricacies of operating HEC-RAS. For those interested in learning more, the HEC-RAS user’s manual is an excellent learning tool and the program is free to the public.

</p><p>The first two figures below are the upstream and downstream water surface profiles modeled by HEC-RAS. There is also a table provided comparing the differences between the profiles estimated by the two different methods at different stations to show consistency between the two methods. While the two different methods modeled similar water surface shapes, the standard step method predicted that the flow would take a greater distance to reach normal depth upstream and downstream of the gate. This stretching is caused by the errors associated with assuming average gradients between two stations of interest during our calculations. Smaller dx values would reduce this error and produce more accurate surface profiles.

</p><p><span typeof="mw:File"></span>

</p><p>The HEC-RAS model calculated that the water backs up to a height of 9.21 meters at the upstream side of the sluice gate, which is the same as the manually calculated value. Normal depth was achieved at approximately 1,700 meters upstream of the gate.

</p><p>HEC-RAS modeled the hydraulic jump to occur 18 meters downstream of the sluice gate.

</p><p><span typeof="mw:File"></span>

</p>

<div class="mw-heading mw-heading2"><h2 id="References">References</h2></div>

/* start https://en.wikipedia.org/ */

.mw-parser-output .reflist{margin-bottom:0.5em;list-style-type:decimal}@media screen{.mw-parser-output .reflist{font-size:90%}}.mw-parser-output .reflist .references{font-size:100%;margin-bottom:0;list-style-type:inherit}.mw-parser-output .reflist-columns-2{column-width:30em}.mw-parser-output .reflist-columns-3{column-width:25em}.mw-parser-output .reflist-columns{margin-top:0.3em}.mw-parser-output .reflist-columns ol{margin-top:0}.mw-parser-output .reflist-columns li{page-break-inside:avoid;break-inside:avoid-column}.mw-parser-output .reflist-upper-alpha{list-style-type:upper-alpha}.mw-parser-output .reflist-upper-roman{list-style-type:upper-roman}.mw-parser-output .reflist-lower-alpha{list-style-type:lower-alpha}.mw-parser-output .reflist-lower-greek{list-style-type:lower-greek}.mw-parser-output .reflist-lower-roman{list-style-type:lower-roman}

/* end https://en.wikipedia.org/ */

</style><div class="reflist">

/* start https://en.wikipedia.org/ */

.mw-parser-output cite.citation{font-style:inherit;word-wrap:break-word}.mw-parser-output .citation q{quotes:"\"""\"""'""'"}.mw-parser-output .citation:target{background-color:rgba(0,127,255,0.133)}.mw-parser-output .id-lock-free.id-lock-free a{background:url("./mw/Lock-green.svg")right 0.1em center/9px no-repeat}.mw-parser-output .id-lock-limited.id-lock-limited a,.mw-parser-output .id-lock-registration.id-lock-registration a{background:url("./mw/Lock-gray-alt-2.svg")right 0.1em center/9px no-repeat}.mw-parser-output .id-lock-subscription.id-lock-subscription a{background:url("./mw/Lock-red-alt-2.svg")right 0.1em center/9px no-repeat}.mw-parser-output .cs1-ws-icon a{background:url("./mw/Wikisource-logo.svg")right 0.1em center/12px no-repeat}body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-free a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-limited a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-registration a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-subscription a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .cs1-ws-icon a{background-size:contain;padding:0 1em 0 0}.mw-parser-output .cs1-code{color:inherit;background:inherit;border:none;padding:inherit}.mw-parser-output .cs1-hidden-error{display:none;color:var(--color-error,#d33)}.mw-parser-output .cs1-visible-error{color:var(--color-error,#d33)}.mw-parser-output .cs1-maint{display:none;color:#085;margin-left:0.3em}.mw-parser-output .cs1-kern-left{padding-left:0.2em}.mw-parser-output .cs1-kern-right{padding-right:0.2em}.mw-parser-output .citation .mw-selflink{font-weight:inherit}@media screen{.mw-parser-output .cs1-format{font-size:95%}html.skin-theme-clientpref-night .mw-parser-output .cs1-maint{color:#18911f}}@media screen and (prefers-color-scheme:dark){html.skin-theme-clientpref-os .mw-parser-output .cs1-maint{color:#18911f}}

/* end https://en.wikipedia.org/ */

</style><cite id="CITEREFUSACE" class="citation web cs1">USACE. "HEC-RAS Version 4.1 User's Manual". Hydrologic Engineering Center, Davis, CA.</cite> <span class="cs1-visible-error citation-comment"><code class="cs1-code">{{cite web}}</code>: </span><span class="cs1-visible-error citation-comment">Missing or empty <code class="cs1-code">|url=</code> (help)</span></span>

</li>

<li id="cite_note-2"><span class="mw-cite-backlink"><b><a href="#cite_ref-2">^</a></b></span> <span class="reference-text"><cite id="CITEREFChaudhry2008" class="citation book cs1">Chaudhry, M.H. (2008). <i>Open-Channel Flow</i>. New York: Springer.</cite></span>

</li>

<li id="cite_note-3"><span class="mw-cite-backlink"><b><a href="#cite_ref-3">^</a></b></span> <span class="reference-text"><cite id="CITEREFMoglen" class="citation web cs1">Moglen, G. <a rel="nofollow" class="external text" href="https://web.archive.org/web/20121105134341/http://filebox.vt.edu/users/moglen/ocf/index.html">"Lecture Notes from CEE 4324/5894: Open Channel Flow, Virginia Tech"</a>. Archived from <a rel="nofollow" class="external text" href="http://filebox.vt.edu/users/moglen/ocf/index.html">the original</a> on November 5, 2012<span class="reference-accessdate">. Retrieved <span class="nowrap">April 24,</span> 2013</span>.</cite></span>

</li>

<li id="cite_note-4"><span class="mw-cite-backlink"><b><a href="#cite_ref-4">^</a></b></span> <span class="reference-text"><cite id="CITEREFChow1959" class="citation book cs1">Chow, V.T. (1959). <i>Open-Channel Hydraulics</i>. New York: McGraw-Hill.</cite></span>

</li>

<li id="cite_note-5"><span class="mw-cite-backlink"><b><a href="#cite_ref-5">^</a></b></span> <span class="reference-text"><cite id="CITEREFChaudhry2008" class="citation book cs1">Chaudhry, M.H. (2008). <i>Open-Channel Flow</i>. New York: Springer.</cite></span>

</li>

<li id="cite_note-6"><span class="mw-cite-backlink"><b><a href="#cite_ref-6">^</a></b></span> <span class="reference-text"><cite id="CITEREFChaudhry2008" class="citation book cs1">Chaudhry, M.H. (2008). <i>Open-Channel Flow</i>. New York: Springer.</cite></span>

</li>

</ol></div></div></div><div><div class="zim-footer">

This article is issued from <a class="external text" title="Last edited on 2024-04-09" href="https://en.wikipedia.org/wiki/?title=Standard_step_method&oldid=1217987743">Wikipedia</a>. The text is available under <a class="external text" href="https://creativecommons.org/licenses/by-sa/4.0/deed.en">Creative Commons Attribution-Share Alike 4.0</a> unless otherwise noted. Additional terms may apply for the media files.

</div>

</div></div>

</div>

</main>

</div>

</body></html>